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Abstract 

We argue that the imaginary parts of the anomalous dimensions 
in the multiparticle sectors of heavy quark effective field theory may 
be removed by a suitable redefinition of the multiparticle states. The 
connection between the imaginary parts of the anomalous dimensions 
and the interquark potential is pointed out. 

1 Introduction 

Heavy quark effective field theory (HqEFT) [p], ||, H, £|, |J has proven to be 
a useful tool for the systematic analysis of systems containing heavy quarks. 
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This approach is based on the heavy mass limit in which the Hilbert space 
decomposes into superselection sectors labelled by the velocities of the heavy 
quarks. Corrections to the leading behavior may be calculated systemati- 
cally and we shall focus in this work on the strong interaction corrections. 
These are treated in renormalization group improved perturbation theory, 
resumming logarithms of the heavy quark masses. Weak interactions will 
introduce operators connecting different velocity sectors and in general their 
anomalous dimensions will depend on the velocities of the heavy quarks. The 
lowest order expressions have been given in || [?]]. 

In the limit itlq — > oo particle and antiparticle number are separately 
conserved and most of the applications discussed up to now deal with the 
one (anti)particle sector of HqEFT. However, in the two particle sector some 
of the anomalous dimensions develop imaginary parts. In the full theory 
imaginary parts of Green's functions emerge from analytic continuation and 
some of these imaginary parts contain logarithms of the masses. In the effec- 
tive theory these logarithms are reproduced by the renormalization group. 
Therefore the anomalous dimensions have to pick up imaginary parts because 
of analyticity of the full theory. 

Complex anomalous dimensions are somewhat unusual and the purpose 
of the present letter is to shed some light on their origin. We shall argue that 
the imaginary parts may be removed by a suitable redefinition of the mul- 
tiparticle states of HqEFT. We shall construct these redefined multiparticle 
states to order a s explicitly in the next section and discuss the extension of 
this method to all orders. 

In the two particle state the imaginary parts of the anomalous dimen- 
sions lead to phases of the Wilson coefficients. Upon closer inspection, these 
phases turn out to be the nonabelian generalizations of the Coulomb phases 
familiar from QED || ||. In a similar way as it may be done for the QED 
case, one may derive an interquark potential suitable for the physics at the 
renormalization scale. This is done up to two loops in section fl[ 

2 Complex Anomalous Dimensions 

Since HqEFT conserves particle and antiparticle number separately, we shall 
begin our discussion in the particle-particle and particle-antiparticle sectors. 
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An effective Hamiltonian 

n 

H eS = Y,Vi(»)Oi(») (i) 

1=1 

having matrix elements with two heavy particles in the final state will be 
expanded in a basis of hermitian operators Oi,i = 1, . . . , n with Wilson co- 
efficients r\i. They are calculated from matching at the heavy quark mass 
scale and subsequent renormalization group running. In such a basis the 
anomalous dimensions will be given by a specific n x n matrix 7. The her- 
miticity of the effective interaction requires real entries for the matrix 7. 
However, in HqEFT imaginary parts appear as coefficients of logarithmic 
divergences as soon as the corresponding operator has matrix elements with 
a two- or more heavy particle state 0. This hints at a possible problem 
for the interpretation of HqEFT for these final states. The imaginary parts 
of these anomalous dimensions have been calculated to leading order in M. 
The generic expression is given by 

w 1 

]mj(a,vv') = -a . =T a ®T a . (2) 

yj(vv') 2 - 1 

The color matrix T a ® T a acts on the two heavy quarks in the final state. 
Coupling the two final state heavy quarks to definite color, it will be sufficient 
to consider the following two sets of operators. In the particle-antiparticle 
sector we have a color singlet and a color octet operator 

J 1 = (h a (v)^h a (v')) (3) 
J 8 < a = {h a {v)^)T:0h p {v')) (4) 

and similarly in the particle-particle sector there is a color antitriplet and a 
color sextet 

J% = ^h a (v))(^(v')) - ^h a (v))(^(v')) (5) 
J 6 a p = $h a (v))($h p (v')) + ($h a (v))($h p (v')). (6) 

In these equations, a, (3 — 1, . . . , 3 and a = 1, ... ,8 are the color indices. The 
Dirac spinors ip and <fi are determined by the specific form of the operators Oi\ 
they may contain fields corresponding to the light degrees of freedom and 
possibly another heavy quark in the initial state. 
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In this basis the color matrix T a ® T a is diagonal which means that the 
imaginary part of the anomalous dimension is diagonal and may be written 

as 




where the color factor Kq is determined by the action of T a ® T a on the 
operators ^,8,3,6; yielding the eigenvalues 4/3, —1/6, 2/3, and —1/3. 

Thus the imaginary parts in fact depend only on the total color of the 
heavy quarks in the final state. This suggests that a suitable redefinition of 
the final state may render the anomalous dimensions real. 



3 Redefinition of the Multiparticle States 

Before we consider the nonabelian case, we shall briefly review the exactly 
soluble abelian case which has been extensively studied in the context of the 
QED infrared problem [10, |{J. The interaction in the abelian case is given 
by 

H I (x )=gJd 3 xj,(x)A^x), (7) 

where the current of a fermion Q is given by 

Ux) = Q(x)^Q{x). (8) 

In the nonrecoil limit corresponding to the heavy mass limit for Q the cur- 
rent (pi) may be rewritten M 



j»(x) -> Jn(x) = jd 3 vv fl 5 3 (x -vx /v )n(v), (9) 
where v = Pq/tuq is the velocity of the heavy particle and 

s)b(v, s) - dt(u, s)d(v, sj) (10) 



n[v) = 

s=± 



is the charge density operator with b(v, s), d(v, s) being the annihilation op- 
erators for particles and antiparticles of spin s moving with velocity v respec- 
tively. The dynamics of the system governed by the Hamiltonian 

Hi(x ) = g J d 3 xJ^x)A»(x) (11) 
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is exactly soluble. The operator U transforming a free state into an interact- 
ing one is given by || 

U = exp(ii2)exp(ifi), (12) 

with the radiation operator 

R = -g J d 4 x&(-x )J fl (x)A^(x) (13) 

and phase operator 

^ = 2" / dAxdA y Q (~ x oM x o - yo) Wx)A»(x) , Ju{y)A"{y)\ , (14) 

This model is exactly soluble because in the abelian case the current (^) 
commutes with itself for all x and y 

[J^x) , J v {y)\ = 0. 

Thus the phase operator may also be written as 

2 

= \ J d A xd A yQ(-x )Q(x - yo)J»(x) [A»{x) , A v {y)\ My). (15) 

Normal ordering of the two currents introduces a mass renormalization term 
which is dealt with as usual. Inserting the c-number commutator of the A 
fields yields for the phase operator 

Q = _ / , = :n(v)n(v'): f (16) 

The divergent xq integration is cut off at small xq by the inverse of the 
mass ttlq of the heavy particle and at large Xq by a scale // which will turn 
out to be the renormalization scale /i. Thus we have 

n = -- / d 3 vd 3 v' , = = :n(?;)n(?/): In f (17) 

2 7 ^(W') 2 - 1 V ^ / 

In the abelian case this operator will contribute to U as a phase factor mul- 
tiplying the multiparticle states. The transformation 
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is unitary by itself and yields the correct Coulomb phases familiar from QED 
if applied to a multiparticle state 0. 

On the other hand, one may also calculate the anomalous dimensions in 
the abelian version of HqEFT with no light quarks which exactly corresponds 
to the above model. From the exact solution (0) we read off that the 
anomalous dimensions are given by the one loop result. The (3 function 
is trivial (/5 = 0) since there is no self coupling of the gauge field and no light 
particles. However, the anomalous dimensions in this simple model develop 
an imaginary part similar to the one in HqEFT which we now trace back to 
the operator Uq. 

In fact, the imaginary part of the anomalous dimensions yields a phase 
factor in the Wilson coefficient. This phase factor for a matrix element M. (vv r ) 
with two heavy particles moving with velocities v and v' respectively in the 
final state is governed by a renormalization group equation with only the 
imaginary part of the anomalous dimension M 

+ ilm^ M(vv') = 0. (18) 

The solution of fll8] ) yields the Coulomb phase which is also obtained from the 
operator Uq acting on the two particle final state. In fact, the operator Uq 
satisfies the renormalization group equation fllBp 

= 0. (19) 

with the imaginary part of the anomalous dimension replaced by the opera- 
tor G 

G = —— f d 3 vd 3 v'^J^^=:n(v)n(v'):. (20) 
2 J J (vv ,)2 _ i y } y } y } 




Note that the hermiticity of G leads to a unitary operator Uq. This makes 
manifest the connection of the Coulomb phases with the imaginary parts 
of the anomalous dimensions. In particular, this argument shows that it is 
appropriate to transfer the phases generated by the imaginary parts of the 
anomalous dimensions from the Wilson coefficients to the states by applying 
the operator Uq to the states. 
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Unfortunately, life is not that simple in the nonabelian case. The currents 
live in the adjoint representation of the gauge group and are given by 

j;(x) = Q a (x)T^ ltl Q p (x), (21) 

where T a is the generator of color SU(3) in the fundamental representation. 
They become in the nonrecoil approximation 



3 fa) - J fa) = Jd d vv,5 d (x-vx /v )n a (v), (22) 

where now 

n a {v) = J2 (bl(v, 8)Tyb p (v, s) - dl(v, s)T«* p d p {v, s)) . (23) 



s=± 



Due to the color structure these currents do not commute any more and an 
exact solution is not possible. Nevertheless we shall proceed along the lines 
suggested by the abelian case and define a nonabelian counterpart of the 
operator G which is given to order a s by 



_ d 3 m 3 v'- r =^=m a (v)n a (v'):. (24) 

2 J J (y V l)2 — i 

This reproduces the imaginary parts of the anomalous dimensions discussed 
in the previous section (cf. (0)) when applied to a two particle state. Fur- 
thermore, Q is a hermitian operator and a renormalization group equation 
analogous to ( pOD may be used to define the nonabelian generalization Uq of 
the operator Uq. However, the self-coupling of the gluons and the existence of 
light quarks imply a nontrivial (3 function. Thus the equation determining 
in the nonabelian case is 

("| +/5 ^ +ie ) Wn = a (25) 

As in the abelian case the operator U is a unitary operator because Q is 
hermitian. The solution of (25) may be obtained in closed form and is given 
by 

U n = exp \- t ^-l f d ^v'^ji==.n\v)n\vyAii(^) \ . (26) 



7 



This operator again defines a unitary transformation on the multiparticle 
states. Accordingly this nonabelian generalization of the Coulomb phase may 
be shifted from the Wilson coefficients to the states, rendering the anomalous 
dimensions real to order a s . 

Unlike the abelian case this is not the end of the story since there are 
higher order contributions to the anomalous dimensions. In general, the 
anomalous dimensions will contain imaginary parts in all orders of a s . Simi- 
lar to the lowest order case one may define an anomalous dimension operator 
by replacing the color factors by appropriate density operators in the follow- 
ing way. The color factor of some set of diagrams is given by a tensor product 
of T a, s corresponding to heavy quark-gluon vertices. However in the nonre- 
coil approximation every T a may be assigned to a fixed velocity v. Thus an 
anomalous dimension operator may be defined by replacing the tensor prod- 
uct of the T a, s by a normal ordered product of the density operators n a (v). 
The generalization of Q to all orders may be obtained by picking out the an- 
tihermitian part of this operator, which exactly corresponds to the imaginary 
part of the anomalous dimension. The hermitian operator Q is now obtained 
by multiplying the antihermitian part of this anomalous dimension operator 
by {-%). 

The kernel Q is put into the renormalization group equation ( p5|) and the 
solution defines a unitary operator Uq which may be used to redefine the 
states, thereby removing the imaginary parts of the anomalous dimension to 
any desired order. 

Q may be expanded in powers of n a (v ) corresponding to two-, three-, . . ., 
n-particle operators 

6 = ^/ d3 ^ d3 ^/(2)K^2):n>i)n 6 (^): (27) 
+ jj d z v x d z v 2 d z v^{v Xl v 2 , v 3 y.n a ( Vl )n b (v 2 )n c (v 3 y. + ■■■. 
The kernels f^" an of these operators are power series in a s with the property 

ygp = 0(a^ 1] ). (28) 
Furthermore, the /fV" 8 ™ may be decomposed into rank n invariant SU(3) 
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tensors; thus the lowest term must be proportional to 5 ab 



vv' 



f$> =-**<*. , ( +•••> ( 29 ) 

J(vv') 2 — 1 

where the dots indicate higher order terms. Up to now only the leading and 
next-to-leading contributions to have been calculated 0, [TTf| . 



4 Potentials and Phases 

Finally we shall point out the relation between the phase operator U^i and 
the nonrelativistic interquark potential. In order to do this we consider a 
final state consisting of a heavy quark and a heavy antiquark in the singlet 
state with velocities v ~ v' . The generalization to other two heavy quark 
final states is obvious since they differ only by color factors. 

We shall consider the nonrelativistic limit, in which vq,v' ~ 1 and we 
have 

' ' (30) 



yjiy ■ v') 2 - 1 



\u\ u 



where u is the velocity of one quark in the rest frame of the other. 

In this limit one obtains a simple result for the imaginary part of the 
anomalous dimension up to two loops for this final state which reads |TT| 



Im7(u) = « s Cf 



a s /31 10 

1 + ^U^-¥ n ' 



where the dots represent terms which are finite asn^O. Cp = 4/3 and Ca = 
3 are the eigenvalues of the Casimir operator in the fundamental and the 
adjoint representation respectively. This imaginary part of the anomalous 
dimension will create a phase of the Wilson coefficient which may be absorbed 
into the states by applying the operator Uq defined in the last section. 

The connection between the potential and the phases appearing in the 
Wilson coefficient proceeds along the lines well known from the QED case P, 
||. Here the Coulomb phases may be traced back to the long range part of 
the potential. In order to connect the phases with the potential we consider 
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a nonrelativistic system of two heavy particles interacting by a spherically 
symmetric potential. The Hamiltonian is 



Separating the cms motion, the time evolution for the relative motion is given 
in the interaction picture by 

U{t, t ) = exp (-i £ dr l/(r 12 (r))) , (33) 

where ^(t) is the operator for the relative coordinate, again in the interac- 
tion picture. 

Since the two particles are very heavy, we can have the particles simulta- 
neously in a state of definite velocity and definite position. We choose 

r 12 (t) = ut, (34) 

where u is the velocity operator. When acting on states containing particles 
with definite velocities, the time evolution becomes a pure phase 

U(t, t ) = exp (-i dr V{ur)^j . (35) 

Introducing the Fourier transform of the potential 

1 f°° 

^( r ) = — T~ I dqqV(q) sin(qr) (36) 
2ir 2 r Jo 

we can rewrite Q55D as 

Cj rt roc qV(q) \ 
-— - / dr dq sin(gwr) . (37) 
Z7T U J to J T J 

For a potential falling off only like l/ri2, the r integration is divergent for t — > 
oo as well as for t — >• 0. However, as argued before, t has to be cut off at 
times of the order of 1/ttiq, since the effective theory is not valid for shorter 
times/larger scales. The other cut off is related to the renormalization point \i 
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by t ~ These cut offs may be shifted from the time integration to the 
integration over the variable q 

[ dr ( dq ^ sm(qur) \— > / dq [ dr ^ sm(qur) (38) 

J t JO T Jl/t=fi JO T 

and the r integration may be performed to yield 

l/(/x; mg) = exp Jj- jT ^ g 2 . (39) 

The exponent is a dimensionless quantity and hence q 2 V(q) is also di- 
mensionless. Calculated from QCD it is thus some function 

<P(a s )=q 2 V(q). (40) 

However, a s is scale dependent and this dependence will render the poten- 
tial nontrivial. Differentiation with respect to the scale /i yields the renor- 
malization group equation (p5|) with the anomalous dimension being purely 
imaginary 

Im 7 (a s ) = -—(f>(a a ). (41) 

This equation finally relates the imaginary part of the anomalous dimension 
to the long range interquark potential for heavy nonrelativistic quarks. 

Including the scale-dependence of oz s , the potential may then be defined 

as 

V(q) = = - lim 4 ™ Im7 ^ (g2)) (42) 

v ' q 2 u^o q 2 v ' 

and we obtain from perturbation theory up to second order in the particle- 
antiparticle singlet channel 



V{q) = 2 



a s (q 2 ) (31 10 

— C A - —n f 



4vr V 9 9 



(43) 



where the two loop running coupling constant has to be used. The QCD 
potential between two static color sources has been calculated long ago JT2], 
|I3]| . Our expression for the potential ([43]), valid in the MS scheme, is in 
complete agreement with these calculations. 
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Finally we point out that for dimensional reasons nontrivial behavior of 
the potential can only be generated by the scale dependence of the strong 
coupling constant. For instance, the linear confinement potential Vn n .(r) oc r, 
which may be obtained from lattice calculations and quarkonia-spectroscopy, 
corresponds to nonperturbative /3-functions. The Fourier-transform of the 
linear potential is 

Vjm.{q) = C~, (44) 
q 2 q 2 

where C is some numerical constant. At this point an additional scale A 2 has 
to enter the game, which has to be generated by dimensional transmutation 
from the scale dependence of the coupling constant. In fact, the general 
renormalization group equation 

d_ 



V 2 — a(/i 2 ) = (3(a(fx 2 )) (45) 



has the solution 




«(^ 2 ) da 



/a(A 2 ) (3(a) J ' 

The potential (04]) can now be expressed in terms of the /3-function 



(46) 



VUq)=C- 2 e W l- -gL) (47) 

q 2 V J*(A2) /3(a) J 



and we can use equation (f42|) backwards to deduce the imaginary part of the 
anomalous dimension corresponding to the linear potential 

C ( f a ^ da \ 
Im7un. = --r^r exp - / — - . (48) 

47T 2 M \ Ja(A 2 ) f3(a) J 

In principle, we could now reconstruct the real part of the anomalous dimen- 
sion from a dispersion relation. However, since the nonrelativistic potential 
approach breaks down for large vv', we do not have any information in this 
region. 



5 Conclusions 

The imaginary parts of the anomalous dimensions of HqEFT may be re- 
moved by a suitable redefinition of the multiparticle states. The redefinition 
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of the multiparticle states is achieved by applying a unitary operator, similar 
to the Coulomb-phase operator in QED. While this phase operator can be 
constructed exactly in the QED case, we can only give a perturbative con- 
struction of the nonabelian analogue of the Coulomb-phase operator. Since 
the phase operator removes the imaginary parts of the anomalous dimen- 
sions, it depends on the renormalization scale; the dependence is governed 
by a renormalization group equation with the imaginary parts of the anoma- 
lous dimensions only. 

Furthermore we have elucidated the connection between the imaginary 
parts of anomalous dimensions and the long range inter-quark potential. The 
potential derived from the HqEFT calculation in the two-particle sector up 
to two loops coincides with earlier perturbative calculations. 

We believe that we have settled the problem of the apparently complex 
anomalous dimensions in the multiparticle sector of HqEFT. The imaginary 
parts are an artifact of an improper definition of the states; after a suitable 
redefinition of the states there is no such thing as a complex anomalous 
dimension. 
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